Draft version February 5, 2008 

Preprint typeset using LSTgX style emulateapj v. 10/09/06 



O 
O 

(N 



INVESTIGATING BINARY PROPERTIES WITH NEXT-GENERATION MICROLENSING SURVEYS 

Cheongho Han 

Program of Brain Korea 21, Institute for Basic Science Research, Department of Physics, 
Chungbuk National University, Chongju 361-763, Korea; cheongho@astroph.chungbuk.ac.kr 

Draft version February 5, 2008 

ABSTRACT 

We explore the usefulness of future gravitational microlensing surveys in the study of binary properties such 
as the binary fraction and the distributions of binary separation and mass ratio by using the binary sample 
detectable through a channel of repeating events. For this, we estimate the rate of repeating microlensing 
eventstoward the Galactic bulge field based on standard models of dynamical and physical distributions of 
Galactic matter combined with models of binary separation and mass function. From this, we find that the total 
number of repeating events expected to be detected from ~ 4-year space-based surveys will be ~ 200—400, 
that is ~ 40-50 times higher than the rate of current surveys. We find that the high detection rate is due to 
the greatly improved sensitivity to events associated with faint source stars and low-magnification events. We 
find that the separation range of the binaries to be covered by the repeating events will extend up to 100 AU. 
Therefore, the future lensing surveys will provide a homogeneous sample that will allow to investigate the 
statistical properties of Galactic binaries unbiased by brightness of the binary components. 
Subject headings: binaries: general - gravitational lensing 
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1. INTRODUCTION 

Massive searches for gravitational microlensing events 
have been and are being carried out by monitoring a large 
number of stars located in the Ga l actic bulge and Mag- 
ellanic Clou ds (lAlcock. et all 11993b lAubourg. et aTJ 120001: 
lUdalskiet al.lll993l) . Current lensing surveys are detecting 
events with a rate of ~ 500 events per year, most of them 
toward the Galactic bulge field, and the total number of de- 
tected events now exceeds 2000. Among them, a significant 
fraction of the events are caused by binaries. 

Although a sizable number of binary lensing events have 
been detected, there have been few investigations about the 
properties of binaries like the binary fraction and the distri- 
butions of the binary separation and mass ratio based on the 
binary event sample. There are several reasons for this. First, 
current lensing detection is strongly biased for binaries with 
separations equivalent to the Einstein radius of the combined 
mass of the binary, and thus the range of binary separation 
covered by lensing is very narrow. Second, for events de- 
tected through this resonant channel, it is difficult to quantify 
the detection efficiency due to the non-linear lensing behav- 
ior. Third, interpreting the lightcurve of these events suffers 
from degeneracy of parameters and thus accurate d etermina- 
tions of the separation and mass ratio are difficult dDominikl 
1999). As a result, there has been no meaningful result about 
the statistical properties of binaries altho ugh there have been 
several reports of binary event samples (|A.lcock. et alj|2000b 
iJarosz vriski 200^ Ijaroszvnski. et alj|2004T) ! 

Besides the resonant channel, binaries can be de- 
tected through a n other channel of repeating events 
dDi Stefano & Mad 1 1996b . Binary detection through 
this channel is possible if an event is produced by a widely 
separated binary lens and the source trajectory approaches 
both lens components. The resulting lightcurve is approxi- 
mated by the superposition of the single lensing events caused 
by the individual binary components. Using the sample of 
binaries detected through this channel has various advantages 
in the studies of binary properties. First, being able to detect 
wide separation binaries, the separation range of binaries that 



can be studied with microlensing can be greatly expanded. 
Second, the detection efficiency is estimated in a straight- 
forward manner and thus can be easily quantified. Third, 
determinations of the binary lens parameters do not suffer 
from degeneracy and thus the parameters can be accurately 
determined. 

Despite these advantages, current lensing surveys are 
very inefficient in detecting binary events through the re- 
peating channel. This is because repeating events require 
special lens-source geom etry and thus they are very rare. 
Di Stefano & Mao ( 1996) pointed out that the detection rate 
of repeating binary lensing events could be increased by con- 
ducting follow-up observations of events. However, these 
follow-ups require long-term observations of the individual 
events. Considering that there are limited resources of follow- 
up observations and the main purpose of the follow-up ob- 
servations cu rrently being con ducted is searchi ng for extra- 
solar planets dBond et al.ll2002l:lPark et al.ll2004l:ICassan et all 
2004), implementing the proposed follow-up observations for 
the detections of wide binaries is practically very difficult. 

Recently, serious discussions are going on about future 
lensing surveys with the prime goal of searching for a large 
sample of extrasolar planets including Earth-mass planets. 
The Microlensing Planet Finder (MPF), which succeeded the 
original concept of the Galacti c Exoplanet Survey Telescope 
(GEST) dBennett & Rhiei r2002). is a space mission exclusive 
for microlensing and it will be equipped with a 1 . 1 m aperture 
telescope, which images a 1.3 deg 2 field of view. The 'Earth- 
Hunter' project is a ground-based survey that plans to achieve 
~ 10 minute sampling by using a network of three 2 m class 
wide field (~ 4 deg 2 ) telescopes scattered over the southern 
hemisphere (A. Gould, private communication). These next- 
generation surveys will have the capability of continuously 
monitoring stars at high cadence by using very large format 
imaging cameras and thus dispense with the need of follow- 
up observations. 

In this paper, we point out the usefulness of the future lens- 
ing surveys in the study of binary properties. We demon- 
strate that a large number of binary events are expected to 
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be detected in the future lensing surveys through the chan- 
nel of repeating events thanks to the greatly improved sen- 
sitivity to events associated with faint source stars and low- 
magnification events. 

The paper is organized as follows. In § 2, we briefly de- 
scribe the basics of repeating binary lensing events. In § 3, 
we estimate the detection rate of repeating binary events ex- 
pected from future lensing surveys based on standard models 
of the physical and dynamical distributions of Galactic matter 
combined with models of binary separation and mass function 
under realistic observational conditions. We also investigate 
the characteristics of the events. In § 4, we discuss about the 
advantages of using the repeating binary event sample in the 
studies of binary properties. In § 5, we summarize the results 
and conclude. 

2. REPEATING EVENTS 

If a lensing event is caused by a binary with a projected 
separation between the lens components significantly larger 
than the Einstein radius of the combined mass of the binary, 
the individual binary components behave as if they are inde- 
pendent single lenses. Then, the resulting lightcurve is well 
represented by the superposition of the single lensing events, 
i.e., 

A~Ai+A 2 -l, (1) 

where A\ and A 2 are the magnifications of the single lens- 
ing events caused by the individual binary components. We 
use the subscripts '1' and '2' to designate the events caused 
by the individual binary components according to the time of 
the source's approach and refer the individual single lensing 
events as the first and second events, respectively. The Ein- 
stein radius is related to the mass of the lens, M, and distance 
to the lens, and source, Ds, by 
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For each single lens event, the magnification is related to the 
lensing parameters by 
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where u is the lens-source separation normalized by the Ein- 
stein radius of each lens component, fg is the timescale re- 
quired for the source to transit the Einstein ring radius (Ein- 
stein timescale), to is the time of the source star's closest ap- 
proaches to the lens, and uq is the lens-source separations at 
that moment (impact parameters) normalized by the Einstein 
radius. The Einstein radius of each binary component is re- 
lated to that of the combined mass by 9ej = y^^E, where 
mi = Mij(M\ +M2) is the mass fraction of the binary compo- 
nent. 

The condition for a wide separation binary to cause a re- 
peating event is that the source trajectory closely approaches 
both lens components. If d is the angular separation between 
the lens components and b\ and bt represent the impact pa- 
rameters (also in angular units) of the source trajectory to the 
individual lens components, then this condition is represented 
by 

^■ rfc " <sin a < '" M+V , (4) 



where a denotes the orientation angle of the source trajectory 
with respect to the binary axis and b t h,i and bth.2 are the max- 
imum values of the impact parameters for the detection of the 
single-lensing events involved with the individual binary com- 
ponents. By normalizing the binary separation in terms of the 
Einstein radius of the combined mass of the binary, s = d/0^, 
and the impact parameters in terms of the Einstein radii of the 
individual lens components, wo.th,; = £>th,//0E,i, equation (01 is 
expressed as 



'«l«0,th.l- V m 2 M 0,th.2 . . „ V m l M 0,th.l + \A«2M0,th.2 

■ < sin a < - — - — . 
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3. DETECTION RATES OF REPEATING EVENTS 



In this section, we estimate the rate of repeating binary lens- 
ing events expected to be detected by future lensing surveys 
and compare the rate to that of current surveys. For this, we 
conduct simulations of repeating events. 

3.1. Lensing Event Simulation 

With the standard models of dynamical and physical distri- 
butions of Galactic matter combined with models of binary 
separation and mass function, we produce repeating events 
toward the Galactic bulge field as follows. 

First, survey is assumed to be conducted toward the Galac- 
tic bulge field and the locations of the source a nd lens are al- 
locate d based on the mass distribution model of H an & Gouldl 
(2003). In the model, the physical distribution of bulge mat- 
ter is based on the deprojected infrared light density profile of 
iDwek et all (1 19951) . specifically model G2 with /? max = 5 kpc 
from their Table 2. The disk matter distribution is modeled by 
a double-exponential disk, which is expressed as 



p(R,z) = />oexp 



r-Ro , \z\ 
h R h z 



(6) 



where (R,z) is the Galactocentric cylindrical coordinates, Ro = 
8 kpc is the distance of the Sun from the Galactic center, po = 
0.06 Mq pc" 3 is the mass density in the solar neighborhood, 
and Iir = 3.5 kpc and h, = 325 pc are the radial and vertical 
scale heights, respectively. 

Second, the velocity distribution of the bulge is deduced 
from the tensor virial theorem, while the disk velocity dis- 
tribution is assumed to have a flat rotation speed of v c = 
220 km s" 1 and velocity dispersions along and normal to the 
disk plane of an = 30 km s" 1 and a± = 20 km s , respec- 
tively. The resulting distribution of the lens-source transverse 
velocity from the combina tions of the bulge an d disk veloc- 
ities is listed in Table 1 of Han & Gould| d 19951) . specifically 
non-rotating barred bulge model. 

Third, the absolute brightnesses of the Galactic bulge 
source stars are assigned on th e basis of the luminosity func- 
tion of Holtzman et al. (1998|), that was constructed from 
Hubble Space Telescope observations. Once the absolute 
magnitude is assigned, the apparent magnitude is determined 
considering the distance to the source star and extinction. The 
extinction is determined such that the source star flux de- 
creases exponentially with the increase of the dust column 
density. The dust column density is computed on the basis of 
an exponential dust distribution model with a scale height of 
h z = 120 pc, i.e. oc exp(-\ z\/h 7 ). Based on the ex tinction deter- 
mination of A v = 1 .28 bv lHoltzman et al.1 dl998l) and the mean 
ratio of extinction to color e xcess of (Rvi) = (Av/E(V— /)) = 
0.964 determined bv lSumil d2004h . we normalize the amount 



HAN 



3 



of /-band extinction so that A/ = 0.49 toward the Baade's Win- 
dow. 

Fourth, the ma ss of the lens i s assigned based on the model 
mass function of iGouldl d2000l) . which is composed of stars, 
brown dwarfs (BDs), and stellar remnants of white dwarfs 
(WDs), neutron stars (NSs), and black holes (BHs). The 
model is constructed under the assumption that bulge stars 
formed initially according to a double power-law distribution 
of the form 



dM 



0.7 M Q 



where 7 = -2.0 for M > 0.7 M Q and 7 = -1 .3 for M < 0.7 M . 
Based on this initial mass function, remnants are modeled by 
assuming that the stars with initial masses 1 M© < M < 8 Mq, 
8 Mq < M < 40 M Q , and m > 40 M Q have evolved into 
WDs (with a mean mass (M) — 0.6 Mq), NSs (with (M) ~ 
1.35 Mq), and BHs (with (M) ~ 5 Mq), respectively. Then, 
the resulting mass fractions of the individual lens components 
are stars : BD : WD : NS : BH = 62 : 7 : 22 : 6 : 3. 

Once the first event is p r oduce d, we then introduce a sec- 
ond event. Following AbJ (1 1 983b . we model the binary sep- 
aration is uniformly distributed in log s, i.e. /(logs) = const. 
To ensure that the part of the lightcurve induced by each of 
the individual lenses can be well approximated by a single 
lens lightcurve, the binary separation is restricted to be larger 
than three times of the combined Einstein ring radius, i.e. 
s > 3. For the mass function of the binary, we test three mod- 
els. In the first model, we assume that the two masses of the 
binary components are drawn independently from the same 
mass function as that of single stars. This model is natural re- 
sults of binary formation where binaries are formed through 
interactions between protostellar disks dPringld 1989b or some 
other forms of capture. We refer this model as the 'capture' 
model. Other possible mechanisms of binary formation are 
fission of a single star and fragmentation of a collapsing ob- 
jects. Numerical calculations suggest that the former process 
results in a mass ratio distribution peaking at around q = 0.2 
( Lucvlll977|). while the l atter results in more equal masses 
(Norman & Wilsonl [l978D. We refer these models as the 'fis- 
sion' and 'fragmentation' models, respectively. In the fission 
model, we model the mass ratio distribution as 



/(<?) = 



f 5^ for q < 0.2 

1-1.25^+1.25 for q> 0.2, 



(8) 



while the distribution of the fragmentation model is modeled 

as 

f(q) = q. (9) 

In the fission model, we set the total mass of the binary to be 
the same as that of a single star because a single star is split 
into two components of a binary in this model. Then, the total 
mass of the binary is drawn from the mass function of single 
stars and the masses of the individual binary components are 
deduced from the mass ratio. In the fragmentation model, on 
the other hand, the mass of the primary is drawn from the 
single star mass function and the mass of the companion is 
determined from the mass ratio. As a result, the average mass 
of the binaries in this model is heavier than that of single stars. 

3.2. Observational Condition and Detection Criteria 

Event rate is not an absolute quantity but is subject to the 
conditions of observation and detection criteria. We estimate 
the rate of repeating events by assuming the following obser- 
vational conditions of the current and future lensing surveys. 



TABLE 1 
Observational Conditions 



specification 


current survey 


future survey 


platform 


ground based 


space mission 


center of field 


(2.0°, -2.7°) 


(1.2°, -2.4°) 


field of view 


21 deg 2 


2.6 deg 2 


photometric precision 


5% at / = 19 


1% at 7 = 21.5 


sampling frequency 


1 time/night 


96 times/night 


fraction of clear nights 


43% 


100% 


observation duration per year 


7 months 


9 months 



NOTE. — Summary of the observational conditions of the cur- 
rent and next-generation microlensing surveys. 

We note that these conditions are based on those of the OGLE 
and MPF surveys, respectively. 
The details of the observational conditions are as follows. 

1 . For the observation fields, we assume that the current 
lensing survey is monitoring toward a bulge field of 
(l,b) ~ (2.0°, -2.7°) with a field of view of (A/, A/3) ~ 
(3.0°, 7.0°). The future survey will monitor toward a 
denser field of (l,b) ~ (1 .2°, -2.4°) but with a narrower 
field of view of (A/, Ab) ~ (0.93°, 2.8°). 

2. For the sampling frequency, we assume that the current 
survey is monitoring the same field with a sampling in- 
terval of 1 day, but the observation is interrupted by spo- 
radic bad weather. The nights of bad weather comprise 
~ 43% of a total 7-month period of the bulge season. 
The future lensing survey will continuously observe the 
target field with a sampling interval of 15 minutes ex- 
cept a 3-month period to avoid the Sun. 

3. Despite that a star is observed by telescopes with sim- 
ilar apertures, the photometric precisions from ground 
and space-based observations are substantially differ- 
ent. This is because the bulge field is very crowded 
and thus ground-based observations suffer from severe 
blending that contributes to the photometric noise. To 
account for blending, we set the photometric precisions 
as 5 % at /= 1 9 for the current survey and 1% at 1 = 21. 5 
for the future survey. 

We summarize the observational conditions in Table Q] 

We judge the detectability of an event based on the un- 
certainties of the lensing parameters recovered from the 
lightcurves. The uncertainties are estimated by computing 
the curvature matrix of \ 2 surface. For the case of a lensing 
lightcurve, the curvature matrix is defined as 



AT*. 



/ _ \ 9F k dF k 1 

b u = L. dpidp . aV 



(10) 



where /V b s is the number of observations, F(t) =A(t)Fs+F\, 
is the observed flux, F& and F\, are fluxes of the unmagni- 
fied source and blended light, a k is the photometric precision 
of each measurement, and pi = (Fs,F\,,uo,tQ,tE) are the five 
lensing parameters required to fit the lightcurve of a stan- 
dard point-source single-lens event. Then, the uncertainties 
of the individual lensing parameters correspond to the diag- 
onal components of the inverse curvature matrix (covariance 
matrix), i.e. 

Op, = Jc~i, c = b~ l . (11) 



Among the lensing parameters, the Einstein timescale is of 
interest because it is the only lensing parameter related to 
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FIG. 1. — Lightcurve of an example lensing event that barely passes the 
imposed detection criteria. See § 3.2, for details of the criteria. 

the physical quantities of the lens system. The timescale 
is also important because the mass ratio between the bi- 
nary lens components of a repeating event is determined by 
q = (/e.2Ae.i) 2 - With the uncertainty of the timescale a ts , we 
set the criteria of event detection as 



f detection 

I non - detection 



ifo> E /f E <0.15, 
otherwise, 



(12) 



where o> E is determined from the lightcurve measured during 
the time of-3?E < t— to < 3?e- We note that the same detec- 
tion criteria are imposed on both the first and second events 
because determinations of the binary parameters such as the 
mass ratio requires not just simple detections of the compan- 
ion signature but equally well-determined timescales of both 
events. We also note that the criteria imposed on the current 
and future surveys are the same for impartial comparison of 
the event rate. In FigureQ] we present an example lightcurve 
of an event that barely passes the imposed detection criteria. 

3.3. Event Rates 

For each event passed the detection criteria, its contribution 
to the total event rate is computed such that 

r Pi ; oc p(D s )p(D L )DlvX, (13) 

where p(D) is the matter density along the line of sight, the 
factor D\ is included to account for the increase of the num- 
ber of source stars with the increase of Ds, v is the lens-source 
transverse speed, and E represents the lensing cross-section. 
The lensing cross-section is proportional to to the physical 
size of the Einstein radius, i.e. £ oc = Z\#e- Since events 
can be detected even when the source does not enter the Ein- 
stein ring, we set the upper limit of the impact parameter large 
enough to account for the contribution of these events to the 
event rate. Although our estimation of the event rate is based 
on numerical calculation to consider detectability depending 
on the observational conditions, we not e that a general ana- 
lytic e xpression for the rate is found in Di Stefa no & Scalzol 
(1999). 

In Table [2] we summarize the resulting rate of repeating 
binary lensing events expected from the current and next- 
generation microlensing surveys under various binary mass 
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FIG. 2. — Distributions of the source brightness of lensing events to be 
detectable from current and future lensing surveys. The distributions with 
light greyscale are for the total single lensing events expected when there 
is no binary lenses. The distributions with dark greyscale are for repeating 
binary events under the binary mass function model where the masses of 
binary components are drawn independently from the same mass function. 
Here we assume a binary fraction of / b j = 0.5. 



TABLE 2 
Comparison of Event Rates 



binary mass 


event rate 


function model 


current 


future 


capture 


5.2/ bi /yr(l.l/ bi %) 


189/ bi /yr(4.3/ bl %) 


fission 


1.9/ bi /yr(0.4/ bl %) 


91/ bi /yr(2.1/ bi %) 


fragmentation 


3.0/ bi /yr(0.6/ bi %) 


128/ bi /yr (2.9/ bi %) 



NOTE. — Comparison of the rates of repeating binary lensing 
events expected from the current and next-generation microlens- 
ing surveys under various models of binary mass function. The 
values in the parentheses are the fractions of repeating events out 
of the total single lensing events expected when there is no binary 
lens. The rate of total single lensing events is ~ 470 events/yr for 
current surveys and ~ 4500 events/yr for future surveys. 



function models. The values in the parentheses are the frac- 
tions of repeating events out of the total single lensing events 
expected when there is no binary lens. We find that the rate 
of total single lensing events is ~ 470 events/yr for current 
surveys and ~ 4500 events/yr for future surveys. From the ta- 
ble, one finds that the rate of repeating events from the future 
survey is ~ 36-47 times higher than that of the current survey 
depending on the adopted models of the binary mass func- 
tion, reaching ~ 91/b;-189/bi per year. Considering 3.7-year 
lifetime of the MPF mission and assuming a binary fraction 
of fb i ~ 0.5 (iDuquennoy & Mayor! Il99lt iFischer & Marcyl 
1992), the total number of binary events expected to be de- 
tected via the repeating channel would be ~ 169-350 events, 
which is large enough for the statistical investigation of the 
binary properties. 
Then, why is the rate of repeating events much higher in 
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FIG. 3. — Distributions of the impact parameters of events detectable from 
the current (upper panel) and future (lower panel) lensing surveys. 

future lensing surveys. We find two reasons for this. The 
first is the greatly enhanced sensitivity of the future survey to 
events associated with faint source stars. This can be seen in 
Figure |2] where we present the distributions of source bright- 
ness of the second part of the repeating events detectable by 
the current and future surveys under the capture model of the 
binary mass function. From the distributions, one finds that 
while the current survey is sensitive only to stars brighter than 
/ ~ 20, the future survey would be sensitive to stars up to 
/ ~ 26. Lensing observation by the current ground-based sur- 
vey is limited by crowding. On the other hand, the space- 
based MPF mission can resolve stars deeper down the lumi- 
nosity function, where there exist much more numerous stars. 
As a result, the event rate of the future survey would be much 
higher than the rate of the current survey despite the coverage 
of substantially smaller field of view. 

Another factor that contributes to the higher detection 
rate of repeating events is the elevated sensitivity to low- 
magnification events. To demonstrate this, we plot the distri- 
butions of the impact parameters in Figure [3] From the com- 
parison of the distributions, one finds that the average impact 
parameter of the events detectable by the future survey is sub- 
stantially larger than the average value of the events detectable 
by the current survey, implying that the future survey will be 
sensitive to lower-magnification events. With the increase of 
the threshold impact parameter, the detection efficiency of re- 
peating events also increases. 

Then, what are the separation range of binaries to be de- 
tected through the channel of repeating events? Another im- 
portant quantity related to the observational side is the wait 
times between the two separate parts of repeating events, f sep . 
In Figure |U we present the distributions of the binary separa- 
tions and wait times. We find that the separation distribution 
peaks at ~ 15 AU, but it extends considerable separations up 
to 100 AU. Therefore, the repeating channel would provide 
sample of binaries with a wide range of separations. The most 
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FIG. 4. — Distributions of the binary separations and wait times for repeat- 
ing events to be detected by future lensing surveys. 

typical wait time is f sep ~ 200 days and extends up to 1000 
days. 

4. DISCUSSION 

Binaries can be detected and characterized by various other 
methods including photometric, spectroscopic, and imaging 
observations of eclipsing, spectroscopic, and visual binaries, 
respectively. Studies of binaries with microlensing have the 
following advantages over these other methods. First, while 
other methods require long-term observations, especially for 
wide-separation binaries, microlensing method can detect and 
characterize binaries from much shorter-term observations. 
Second, unlike other methods, microlensing method has no 
bias against dark or faint binary components. Third, while 
most other methods can be applicable only to binaries located 
relatively at small distances from the Earth, binaries to be de- 
tected by lensing surveys will be located at Galactic-scale dis- 
tances. Fourth, next-generation lensing surveys will be able to 
produce a homogeneous sample of large binaries. Therefore, 
microlensing technique will be able to provide better statistics 
of Galactic population of binaries. 

5. CONCLUSION 

We explored the usefulness of future lensing surveys in the 
study of binary properties by using the binary sample de- 
tectable through a channel of repeating events. To demon- 
strate this, we estimated the rate of repeating events based on 
standard models of mass function and dynamical and physical 
distributions of Galactic matter combined with models of bi- 
nary separation and companion mass function. From this, we 
found that the total number of repeating events expected to be 
detected from ~ 4-year space-based survey would be ~ 200- 
400, which is ~ 40-50 times higher than the rate of current 
surveys, and the range of binary separations to be covered by 
lensing extends to up to 100 AU. We found that the high de- 
tection rate is due to the greatly improved sensitivity of the 
future survey to events associated with faint source stars and 
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low-magnification events. Therefore, the future lensing sur- 
veys will provide a homogeneous binary sample that will al- 
low to investigate the statistical properties of Galactic binaries 
unbiased by the brightness of the binary components. 



This work was supported by the grant (C00072) of the Ko- 
rea Research Foundation. 
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